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DYNAMICAL SYSTEMS. 
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Using the tools of Differential Geometry, we define a new fast chaoticity indica- 
tor, able to detect dynamical instability of trajectories much more effectively, (i.e., 
quickly) than the usual tools, like Lyapunov Characteristic Numbers (LCN's) or 
Poincare Surface of Section. Moreover, at variance with other fast indicators pro- 
posed in the Literature, it gives informations about the asymptotic behaviour of 
trajectories, though being local in phase-space. Furthermore, it detects the chaotic 
or regular nature of geodesies without any reference to a given perturbation and 
it allows also to discriminate between different regimes (and possibly sources) of 
chaos in distinct regions of phase-space. 

Chaotic dynamics is believed to be the rule rather than the exception in any- 
generic nonlinear dynamical system; however, there is no general fast method able 
to determine what's the fraction of the allowed phase space occupied by chaotic 
orbits. This problem exists in the case of few degrees of freedom (dof ) systems and 
for many dimensional ones as well, though it .originates from different causes in the 
two situations. Recently the authors deriveda a geometric indicator of Chaos, able 
to single out the chaotic or regular character of individual orbits simply computing 
a pair of short time correlation functions of geometric quantities defined on the 
tangent bundle of the system. Here we discuss some points which, in our opinion, 
make our indicator preferable with respect to other "fast" indicators, and, from a 
practical viewpoint, we compare its performances with those of the others, showing 
the faster convergence properties and the richer amount of information it contains. 

In addition to the well known reduction of dynamics to geometry which-.can be 
accomplished on the grounds of the Maupertuis principle of Least ActionEJ, there 
is a variety of possibilities of transcription of dynamical trajectories into geodesic 
flows over suitable manifolds. A particularly appealing choice, in the case of generic 
Lagrangian systems, with a more general than simply quadratic dependence on ve- 
locities, is given by the Finsler approach, which furnishes a metric on the tangent 
bundle of the Lagrangian system under examinationtlQ. Dynamical instability is 
linked to curvature properties of the manifold through Jacobi-Levi-Civita equation 
for geodesic spread (or a suitable generalization thereof), which describes the evo- 
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lution of a generic perturbation Sq a to a given geodesic, as a function of the affinc 
parameter, s, along the geodesic itself: 



V fVSq a 



ds V ds 



«W = o, (l) 



where the stability tensomu, H. a c , defined along any geodesic, is related to the 
Finsler curvature tensor, K a f, c d, associated to the metric, by 

U a c = K a bcd x' b x' d , (2) 

and x' is the unit tangent vector to the geodesic. Finsler curvature tensor (and con- 
sequently the stability tensor) is a pretty complicated object Em; however the main 
features of the geodesic flow depend essentially on suitable averages and moments 
of scalars built from TL a c . In the case of iV-dof systems, the curvature properties 
of the associated .Finsler manifold are described by the N non trivial principal sec- 
tional curvaturesu, {Xi(s)}, which are specific of each given path, i.e., depend not 
only on the point on the configuration space but also on the tangent vector along 
the given geodesicQ- If we restrict to two dof systems, the curvature properties of 
the Finsler manifold are then described by the two principal sectional curvatures 
{Ai, A2}. The mean curvature n and the anisotropy $ are defined, in terms of the 
principal sectional curvatures, simply as 

. X 1 + A 2 _ Tr(H) Ric F (x') X 1 - A 2 

k = = = ; v — . (0) 

2 2 2' 2 V ' 

where it has been exploited the link between the (Finsler generalization of) Ricci 
curvature along the flow and the trace of the stability tensor. For any realistic 
dynamical system, it is easy to verify that the associated Finsler manifold and the 
geodesies on it are characterized by principal sectionals which fluctuate around a 
well defined average^, Ai = \i{s). Schur theoremlMJ establishes a link between the 
anisotropy of a manifold and the fluctuations of the curvatures moving along a given 
path, in particular along a geodesic, and suggests a possible tool to characterize 
the interplay between them. 

Introducing the correlation function of any two given observables, X{s) = 
X[q(s),q'(s)] and Y(s) = Y[q(s), q'(s)] computed over a finite interval, S, along 
a geodesic with initial conditions (q o ,q' ) = [q(0), q'(0)], 

Cs[X,Y] = {X ' Y)s , (4) 

«* 2 ) S • <r 2 )s) 1/2 

we define the functional Rf[S] over the given segment of the geodesic, as 

R F [S] = R F (q ,c{ a ) = J-^-j > 0, (5) 



"This is one of the main differences which make the Finsler geometrization more adeherent to the 
complete dynamics than others, like Jacobi or Eisenhart. 

b The only requirement is that the phase space is bounded and the dynamical system is singularity- 
free; though with some cautions the application of our approach can be extended even to singular 
lagrangians. 
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where 6X(s) = X(s) - (X). We have showD&B that the stability properties of 
the geodesic flow, and then the degree of chaoticity of the trajectories, possess an 
one-to-one correspondence with the value assumed by this functional over any short 
segment of a geodesic. Just by computing Rp over a few periods along an orbit, 
suffices to have a reliable indication on its regular or chaotic character. Here we 
report some supplementary results, still referred to the paradigmatic Henon-Heiles 
systemEj, whose Lagrangian is 



which emphasize the much faster convergence of our indicator with respect to the 
usual ones. A pictorial taste of the reliability and effectiveness of the Finsler in- 
dicator, has been given, for example, in figures 1 and 2 of ref.u, where we made a 
synoptic comparison of the gray-scale map of Rf with the PSS for a typical energy 
value (E = 0.125), for which the phase space of this system is almost equally popu- 
lated by chaotic and regular trajectories; noticing that the integration time S, used 
to produce the map of Rp [S^-map is there about two orders of magnitude shorter 
than the one required to obtain the PSS with a comparable resolution level. Fur- 
thermore, we observe that, if a better resolution of the substructures on the surface 
is sought, the computation of Rf requires obviously only a more careful numerical 
integration, whereas the PSS needs also to improve the bisection algorithm needed 
to find the exact intersection point. More detailed results and discussions, and ex- 
tensive applications to other few dimensional systems can be found elsewhereQ. To 
further evaluate quantitatively the fast convergence of the proposed indicator, we 
report in figure 1 the values of the LCN's and of Rf[S] as a function of £im^] for 
some typical orbits of the same system, which makes self-evident the much faster 
signature of regular or chaotic motion obtained using the geometric indicator. In 
addition, we stress also that the values of our indicator are completely independent 
from any given choice of the initial perturbation vector, whereas this is not true for 
the computation of the LCN's, whose values, especially for few dof systems, can be 
strongly influenced by the choice of the orientation of the initial disturbance vec- 
tor. This is mostly evident if the reference orbit is near the border between regular 
islands and the stochastic sea, where stickness phenomena are expected to occur. 
The three plots in the figure indeed demonstrate all what has been said here: 

• Plot [A], showing the long time behaviour of LCN's, singles out that for some 
rather sticky orbits, the time required for the non vanishing LCN's to converge 
towards the true asymptotic values can be about three orders of magnitude 
larger than the associated Lyapunov time, tl ~ LCN^ 1 . 

• Moreover, the two curves labeled as (3) and (36) show that up to t = 10 3 the 
behaviour of two LCN's, computed along the same orbit but with a different 
choice of the orientation of the perturbation vector in phase space, can display 
completely different behaviours. 

c See iEj for the relationship between the Newtonian time t and the Finsler afBne parameter sp. 
We recall only that we can always choose, as we did here, sp(t) = t. 
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Figure 1. Convergence properties of LCN's and of the (normalized) geometric indicator Rp: panels 
[A] and [B] display respectively the long and short time behaviours of LCN's for a set of initial 
conditions on the energy surface E = 0.125. Panel [C] shows the values of Rp for the same 
trajectories and for two further regular orbits. The much faster convergence of the latter indicator 
is evident. 



• This is confirmed by the behaviour reported in the plot [B], where the same 
curves of the previous graph are restricted to a relatively short interval, which 
reveals how, up to t ~ 200, LCN's are completely inadequate to discriminate 
between regular and chaotic orbits. 

• Plots in [C] show instead that the (normalized!!) geometric indicator Rp possess 
much better convergence properties: it is evident that already at t ~ 50 the 
discrimination between chaotic and regular orbits is clear and the fluctuations 
around the average for regular orbits are almost negligible. Moreover the values 
assumed by our indicator give a measure about the degree of regularity of non 
chaotic trajectories, indicating which topological region of phase space will 
become chaotic earlier, as we increase the perturbation parameter, i.e., in this 
case, the energy. 

A couple of final comments are in order. 
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The first one addresses the dichotomy encountered in the characterization of dy- 
namical systems, where no intermediate status exist between ordered and chaotic 
motions, for a given dynamical system: the picture based on the usual tools shows 
that the change from asymptotically ordered to chaotic dynamics, presents itself 
indeed as a phase transition, along with a variation of an external parameter (e.g., 
the energy), with a sudden jump of asymptotic observables (e.g., LCN's) from a 
vanishing value to a finite one. The proposed geometric indicator, though confirm- 
ing the occurrence of a change of state phenomenology, brought to evidence even by 
local correlation functions, indicates however the presence of gradual, premonitory 
changes of local observables, accompanying the approach to the critical point. 

The second remark deals with some recently pjpposed so-called fast indicators, 
like rotation angles and Fast Lyapunov IndicatorsB'tll. These interesting alternative 
tools to detect chaoticity, when applied to continuous time dynamical systems, 
give "fast" informations about the "local" dynamics; if, inseatd, an asymptotic 
characterization is sought, then their computation needs, in general, as much as 
integration time as that required by the evaluation of LCN's. In addition they also 
make necessarily reference to a given choice of the perturbation vector. Nevertheless 
they are best suited for the study of the KAM and Nekhoroshev thresholds in the 
case of discrete time maps, where our approach is instead unsuitable. 

In any case, we point out that our geometric indicator, which is defined only in 
terms of the region of the tangent bundle explored by a trajectory, gives a deeper 
conceptual insight on the sources of Chaos in dynamical systems, just because is 
strictly linked to the properties of the underlying manifold and make no reference 
to any perturbations. 
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